Modeling coronal magnetic field using spherical geometry: cases with several active regions

The magnetic fields in the solar atmosphere structure the plasma, store free magnetic energy and produce a wide variety of active solar phenomena, like flare and coronal mass ejections(CMEs). The distribution and strength of magnetic fields are routinely measured in the solar surface(photosphere). Therefore, there is considerable interest in accurately modeling the 3D structure of the coronal magnetic field using photospheric vector magnetograms. Knowledge of the 3D structure of magnetic field lines also help us to interpret other coronal observations, e.g., EUV images of the radiating coronal plasma. Nonlinear force-free field (NLFFF) models are thought to be viable tools for those task. Usually those models use Cartesian geometry. However, the spherical nature of the solar surface cannot be neglected when the field of view is large. In this work, we model the coronal magnetic field above multiple active regions using NLFFF extrapolation code using vector magnetograph data from the Synoptic Optical Long-term Investigations of the Sun survey (SOLIS)/ Vector Spectromagnetograph (VSM) as a boundary conditions. We compare projections of the resulting magnetic field lines solutions with their respective coronal EUV-images from the Atmospheric Imaging Assembly (SDO/AIA) observed on October 11, 2011 and November 13, 2012. This study has found that the NLFFF model in spherical geometry reconstructs the magnetic configurations for several active regions which agrees with observations. During October 11, 2011 observation, there are substantial number of trans-equatorial loops carrying electric current.Comment: 3 Figures, Submitted to Astrophysics and Space Science Journa

Electron impact polarization expected in solar EUV lines from flaring chromospheres/transition regions

We have evaluated lower bounds on the degree of impact Extreme Ultraviolet/Ultraviolet (EUV/UV) line polarization expected during solar flares. This polarization arises from collisional excitation by energetic electrons with non-Maxwellian velocity distributions. Linear polarization was observed in the S I 1437 A line by the Ultraviolet Spectrometer and Polarimeter/Solar Maximum Mission (UVSP/SMM) during a flare on 15 July 1980. An early interpretation suggested that impact excitation by electrons propagating through the steep temperature gradient of the flaring transition region/high chromosphere produced this polarization. Our calculations show that the observed polarization in this UV line cannot be due to this effect. We find instead that, in some flare models, the energetic electrons can produce an impact polarization of a few percent in EUV neutral helium lines (i.e., lambda lambda 522, 537, and 584 A)

Intercomparison of numerical models of flaring coronal loops

The proposed Benchmark Problem consists of an infinitesimal magnetic flux tube containing a low-beta plasma. The field strength is assumed to be so large that the plasma can move only along the flux tube, whose shape remains invariant with time (i.e., the fluid motion is essentially one-dimensional). The flux tube cross section is taken to be constant over its entire length. In planar view the flux tube has a semi-circular shape, symmetric about its midpoint s = s sub max and intersecting the chromosphere-corona interface (CCI) perpendicularly at each foot point. The arc length from the loop apex to the CCI is 10,000 km. The flux tube extends an additional 2000 km below the CCI to include the chromosphere, which initially has a uniform temperature of 8000 K. The temperature at the top of the loop was fixed initially at 2 X 1 million K. The plasma is assumed to be a perfect gas (gamma = 5/3), consisting of pure hydrogen which is considered to be fully ionized at all temperatures. For simplicity, moreover, the electron and ion temperatures are taken to be everywhere equal at all times (corresponding to an artificially enhanced electron-ion collisional coupling). While there was more-or-less unanimous agreement as to certain global properties of the system behavior (peak temperature reached, thermal-wave time scales, etc.), no two groups could claim satisfactory accord when a more detailed comparison of solutions was attempted

Simulations of Astrophysical Fluid Instabilities

We present direct numerical simulations of mixing at Rayleigh-Taylor unstable interfaces performed with the FLASH code, developed at the ASCI/Alliances Center for Astrophysical Thermonuclear Flashes at the University of Chicago. We present initial results of single-mode studies in two and three dimensions. Our results indicate that three-dimensional instabilities grow significantly faster than two-dimensional instabilities and that grid resolution can have a significant effect on instability growth rates. We also find that unphysical diffusive mixing occurs at the fluid interface, particularly in poorly resolved simulations.Comment: 3 pages, 1 figure. To appear in the proceedings of the 20th Texas Symposium on Relativistic Astrophysic

Three-dimensional adaptive evolution of gravitational waves in numerical relativity

Adaptive techniques are crucial for successful numerical modeling of gravitational waves from astrophysical sources such as coalescing compact binaries, since the radiation typically has wavelengths much larger than the scale of the sources. We have carried out an important step toward this goal, the evolution of weak gravitational waves using adaptive mesh refinement in the Einstein equations. The 2-level adaptive simulation is compared with unigrid runs at coarse and fine resolution, and is shown to track closely the features of the fine grid run.Comment: REVTeX, 7 pages, including three figures; submitted to Physical Review

The Dynamic Formation of Prominence Condensations

We present simulations of a model for the formation of a prominence condensation in a coronal loop. The key idea behind the model is that the spatial localization of loop heating near the chromosphere leads to a catastrophic cooling in the corona (Antiochos & Klimchuk 1991). Using a new adaptive grid code, we simulate the complete growth of a condensation, and find that after approx. 5,000 s it reaches a quasi-steady state. We show that the size and the growth time of the condensation are in good agreement with data, and discuss the implications of the model for coronal heating and SOHO/TRACE observations.Comment: Astrophysical Journal latex file, 20 pages, 7 b-w figures (gif files